Superconducting "metals" and "insulators" 
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We propose a characterization of zero temperature phases in disordered superconductors on the 
basis of the nature of quasiparticle transport. In three dimensional systems, there are two distinct 
phases in close analogy to the distinction between normal metals and insulators: the superconduct- 
ing "metal" with delocalized quasiparticle excitations and the superconducting "insulator" with 
localized quasiparticles. We describe experimental realizations of either phase, and study their gen- 
eral properties theoretically. We suggest experiments where it should be possible to tune from one 
superconducting phase to the other, thereby probing a novel "metal-insulator" transition inside a 
superconductor. We point out various implications of our results for the phase transitions where 
the superconductor is destroyed at zero temperature to form either a normal metal or a normal 
insulator. 



I. INTRODUCTION 

While Cooper pairing in a superconductor is usually 
associated with a gap in the energy spectrum, it is quite 
well-known that the gap is by no means indispensable. 
Indeed, there are a number of situations in which super- 
conductivity occurs with no gap in the quasiparticle ex- 
citation spectrum. This possibility was first discussed in 
pioneering work by Abrikosov and Gorkov |l| for s- wave 
superconductors in the presence of magnetic impurities. 
Another novel instance is provided by Type // s-wave 
superconductors in strong magnetic fields close to, but 
smaller than, Hc2- Yet another example of gapless su- 
perconductivity, of much current interest, is the d^^^y^ 
superconductor; the d-wave characteristics of the pairing 
is reflected in four nodes in the gap function along which 
gapless quasiparticle excitations result. In all cases, the 
presence of low energy quasiparticle states has a profound 
effect on the low temperature thermodynamic and trans- 
port properties of the superconductor. 

The effect of static disorder introduced by frozen impu- 
rities in such a gapless superconductor raises a number 
of fundamental questions. The low energy quasiparti- 
cles can either be delocalized and free to move through 
the sample, or can be localized by the disorder. These 
two possibilities correspond to two distinct superconduct- 
ing phases that are distinguished by the nature of quasi- 
particle transport. This distinction has a direct analogy 
in the physics of normal (non-superconducting) systems 
where again there are two possible phases distinguished 
by the nature of transport - the metal with diffusive 
transport at the longest length scales, or the insulator 
characterized by the absence of diffusion. Recent work 
P,B| has addressed the possible existence and properties 
of these superconducting phases in the context of dirty 
dx2-y2 superconductors. In particular, it was shown B] 
that quantum interference effects destabilize the super- 
conducting phase with delocalized quasiparticles in two 
dimensions. Instead, the quasiparticle excitations in a 



two dimensional superconductor are generically always 
localized M. In three dimensions however, both kinds of 
superconducting phases are stable. 

It is clear that the issues raised above are germane to 
the properties of all superconducting systems. Consider, 
for instance, the case of dirty Type // s-wave supercon- 
ductors in strong magnetic fields. For fields above Hd, 
and in the absence of impurities, such a superconductor 
exists in the Abrikosov vortex lattice phase. The quasi- 
particle excitation spectrum in the presence of a single 
isolated vortex (relevant when the field is just over Hd) 
was shown by Caroli et. al. [g| to consist, at low ener- 
gies, of a set of discrete energy levels corresponding to 
states bound to the core of the vortex. The lowest lying 
excited state is separated from the ground state by an 
energy gap that is much smaller than the bulk gap, but 
is nevertheless non-vanishing. Impurities alter this den- 
sity of states, and could potentially close the (mini)gap. 
Furthermore, the quasiparticle states are no longer ex- 
tended even along the direction of the vortex line due to 
the strong effects of localization by the disorder in one di- 
mension [^,§,^ . Impurities also destroy the translational 
order of the Abrikosov vortex lattice, pinning the vortices 
into a glassy phase. When many vortices are present, the 
quasiparticles can tunnel from the core of one vortex to 
the other. The amplitude to tunnel depends on the over- 
lap of the wavefunctions of the core states corresponding 
to the individual vortices. This in turn is determined by 
the density and spatial configuration of the vortices. At 
low fields, the tunneling is insignificant, and the quasi- 
particle states are localized. 

With increasing field however, the density of vortices, 
and hence the tunneling, increases. The nature of the 
quasiparticle states - i.e whether they are delocalized or 
localized - is determined by the interplay between tun- 
neling and the disorder. One can envision two scenarios 
for what happens as the field is increased toward Hc2 at 
T = 0. One possibility is that the quasiparticle states 
remain localized all the way upto Hc2- We suggest that 



this is the case if the non-superconducting ground state 
obtained for H > Hc2 is a locahzed insulator. A second 
possibiUty is that the quasiparticle states undergo a delo- 
cahzation transition at some field Hd < Hc2 which there- 
fore precedes the destruction of superconductivity. We 
suggest that this happens if the ground state at iJ > Hc2 
is a normal metal. 
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FIG. 1. A possible phase diagram for a dirty Type II su- 
perconductor 



Similar considerations apply as well to the situation 
where the superconductivity is destroyed at zero tem- 
perature and zero magnetic field, say, by increasing the 
amount of impurity scattering. If the transition is to a 
normal metal, it is likely to be from a superconducting 
phase with delocalized quasiparticle excitations. On the 
other hand, far from the transition in the superconduct- 
ing phase, the quasiparticle states at the Fermi energy, if 
any, will be strongly localized. We expect therefore that a 
delocalization phase transition from this phase to the su- 
perconductor with extended quasiparticle states precedes 
the ultimate transition from the superconductor to the 
normal metal. Again, if the destruction of superconduc- 
tivity leads to an insulating phase, it is natural to expect 
that the transition is directly from the superconductor 
with localized quasiparticle excitations. 

Thus the two kinds of superconducting phases, distin- 
guished by the nature of quasiparticle transport, are both 
realizable in a variety of situations. Despite the analogies 
with normal metals and insulators, as pointed out in Ref. 
H the problem of quasiparticle transport in a dirty su- 
perconductor is conceptually very diff'erent from the cor- 
responding problem in a normal metal. This is because, 
unlike in a normal metal, the charge of the quasiparticles 
in a superconductor is not conserved, and hence cannot 
be transported through diffusion. However, the energy 
of the quasiparticles is conserved, and in principle the 
energy density could diffuse. In addition , in a singlet 
superconductor (in the absence of spin-orbit scattering) 
the spin of the quasiparticles is a good quantum num- 



ber and is conserved. Thus the quasiparticle spin den- 
sity could also diffuse in the superconductor. We may 
therefore characterize the two distinct kinds of supercon- 
ductors by means of their energy (and spin) transport 
properties. The phase with delocalized quasiparticle ex- 
citations has diffusive transport of energy (though not of 
charge) at the longest length scales. It is hence appro- 
priate to call it a "superconducting thermal metal" . In 
cases where the quasiparticle spin is conserved, this phase 
also has spin diffusion at the longest length scales. On 
the other hand, in the phase with localized quasiparticle 
excitations, there is no diffusion of energy or spin at the 
longest length scales. Hence such a phase may be called 
a "superconducting thermal insulator" . 

For the sake of concreteness, we focus for most of the 
paper on the properties of dirty Type II superconductors 
in strong magnetic fields. We discuss the possible phase 
diagrams, and analyse the properties of the two super- 
conducting phases. This is first done in a BCS model 
of non- interacting quasiparticles, and assuming spin ro- 
tational invariance. In striking contrast to normal (non- 
superconducting) systems, quantum interference effects 
lead to singular corrections to the quasiparticle density 
of states in the superconductor Q . The density of states 
at the Fermi energy is non-zero in the superconducting 
"metal" phase, but has a ^/E cusp as a function of energy 
E (measured from the Fermi energy). The superconduct- 
ing "insulator" phase, on the other hand, has a vanishing 
density of states at the Fermi energy. We will provide 
supporting numerical evidence for this result in agree- 
ment with earlier analytical calculations [H . Interactions 
between the quasiparticles could potentially modify the 
properties of either phase in important ways. These are 
analysed next. In the "metal" phase, we show the exis- 
tence of effects analogous to the Altshuler-Aronov singu- 
larities in a normal metal due to the interplay of diffusion 
and interactions. Interaction effects are more crucial in 
the "insulator" - we provide simple arguments showing 
that arbitrarily weak repulsive interactions lead to the 
formation of free paramagnetic moments. The ultimate 
fate of the "insulator" in the presence of these free mo- 
ments is a complicated problem, and we will not address 
it here. Weak attractive interactions however appear in- 
nocuous. We then conclude by discussing the many im- 
plications of this paper for experiments on various super- 
conducting systems. 



II. MODELS 



In this section, we will introduce models appropriate 
to describe the physics of the quasiparticles in a dirty su- 
perconductor in the pinned vortex phase. We begin with 
a general BCS quasiparticle Hamiltonian: 
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Here m is the mass of the quasiparticles, and Ep is the 
Fermi energy. The function V{x) is a random potential 
due to impurities that scatter the quasiparticles. The 
physical magnetic field is introduced through the vector 
potential A. In the pinned vortex phase of a dirty super- 
conductor, the gap function A{x) may be considered a 
random complex function of the position x due to the ran- 
dom positions of the vortices. Both the vector potential 
and the complex gap function break time-reversal sym- 
metry. As the phase of A(x) winds by 2tt on encircling a 
vortex, the corresponding term in the Hamiltonian that 
breaks time reversal symmetry is of order the zero field 
gap Aq. 

The spinful quasiparticles experience a shift in energy 
due to Zeeman splitting (not included in the Hamilto- 
nian Eqn.|l|). The largest field that we consider, Hc2, is 
of order — /Q^ where ^o is the coherence length, making 
the associated Zeeman energy E^m of order /j^/m^^. As 
fo is related to Aq by Aq 



^, we have 
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Here pp = mvp is the Fermi momentum. The last esti- 
mate uses the fact that the coherence length is typically 
a few thousand times larger than the Fermi wavelength 
in conventional low- Tc superconductors. Thus the Zee- 
man energy is negligible compared to the (zero field) gap, 
and we will ignore it for most of our discussion. We 
will also, for the most part, assume the absence of any 
strong spin-orbit scattering. The system is then spin ro- 
tational invariant. The Hamiltonian above also describes 
non-interacting quasiparticles. Interaction effects are po- 
tentially quite important, and will be discussed at some 
length later on in the paper. 

It will often be convenient to think in terms of a lattice 
version of the Hamiltonian Ec^n. n. 
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The Hamiltonian, Eq. ||, then takes the form 



A* _.* ) dj — 2_^"'iHijdj. (5) 



Note that SU{2) spin rotational invariance requires 



cTyH,jay - -iJ*. 



(6) 



The advantage of going to the d representation is that 
the Hamiltonian conserves the number of d particles. The 
transformation Eq. implies that the number of d par- 
ticles is essentially the z component of the physical spin 
density: 
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A spin rotation about the z axis corresponds to a C/(l) 
rotation of the d operators. This U{1) is clearly present 
in the d Hamiltonian. Invariance under spin rotations 
about the x or y axes is however not manifest. 

The Hamiltonian may be diagonalized by the Bogoli- 
ubov transformation 



n 



(8) 
(9) 



Here the "fn],Jni are fermionic operators, and the func- 
tions Un(i),Vn{i) are determined by solving the eigen- 
value equation: 
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Note that these are just the familiar Bogoliubov- 
deGennes equations [l|. The symmetry Eqn. ra which 
follows from physical spin rotation invariance then im- 
plies that 
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(11) 



Hermiticity implies the condition tij = t*^, and spin ro- 
tation invariance requires Ay = Aj^. Note that this is 
the most general spin rotational invariant Hamiltonian 
describing a superconductor with broken time reversal 
symmetry. 

For some purposes, it is useful to use an alternate rep- 
resentation in terms of a new set of d-operators defined 
by: 



H] — ^i] ; 
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is an eigenfunction of H with eigenvalue — £"«. Thus the 
eigenvalues come in pairs (£"«, —En). In terms of the 7„ 
operators, the Hamiltonian becomes 



^L = E -^" (7icr^7„ 



(12) 



The ground state of the superconductor is obtained by 
filling all the negative energy levels, i.e by occupying all 



(4) the states created by 7, 



"i- 



A. Phase diagrams 

In this subsection, we discuss the zero temperature 
phase diagram of the dirty Type // superconductor in the 
mixed state when the vortices are pinned by the disorder. 
We are interested in characterizing the nature of quasi- 
particle transport in such a superconductor. By anal- 
ogy with normal (non-superconducting) systems, we ex- 
pect to distinguish two qualitatively different situations. 
The quasiparticle wavefunctions (m„, Vn) for states at the 
Fermi energy (if any) may either be extended through 
out the system, or they may be localized. Again, in anal- 
ogy with normal metals, we expect that the phase with 
extended quasiparticle wavefunctions is characterized by 
diffusive transport. However, as emphasized in Ref. 0, 
the electric charge of the quasiparticles is not conserved 
by the BCS Hamiltonian, and hence cannot diffuse. The 
energy and spin densities of the quasiparticles are still 
conserved quantities and are thus capable of being trans- 
ported through diffusion. Thus the superconductor with 
extended quasiparticle states at the Fermi energy is char- 
acterized by diffusive transport of energy and spin at the 
largest length scales. We will refer to this ground state as 
a "superconducting thermal metal" . Similarly, the phase 
with localized quasiparticle wavefunctions at the Fermi 
energy is characterized at T = by the absence of diffu- 
sion of energy and spin . We will refer to this phase as 
the "superconducting thermal insulator" . 

Recent studies of Anderson localization in a supercon- 
ductor with spin rotation invariance show that J2l in two 
or lower dimensions, the quasiparticle wavefunctions are 
generically always localized. Above two dimensions how- 
ever, phases with extended or localized states are possi- 
ble. While these results are similar to the corresponding 
results in a normal metal, there is a striking difference 
|g|: the quantum interference effects leading to localiza- 
tion give rise to singular corrections in the quasiparticle 
density of states unlike in a normal metal. We will pro- 
vide additional support for this result, and explore its 
consequences in later sections. 

First consider the situation at low fields H just bigger 
than Hci when the inter-vortex spacing is large. For a 
single isolated vortex, in the absence of impurities, there 
exist states bound to the core B whose energies form a 
discrete set of minibands. (The (mini)bandwidth is en- 
tirely due to motion along the vortex line). There still is 
a gap to all excitations of order -g^. This is very small 
compared to Aq, which in turn is of order ksTc- In the 
presence of impurities, the quasiparticle wavefunctions 
are localized along the direction of the vortex line P,0,p| • 
Further, for large inter- vortex spacing , quasiparticle tun- 
neling from the core of one vortex to another is insignifi- 
cant. Thus we expect no energy or spin transport by the 
quasiparticles at the longest length scales; the system is 
in the thermal insulator phase. 



Upon increasing the magnetic field towards Hc2, it 
becomes important to include tunneling from the core 
of one vortex to the other. If the tunneling is suffi- 
ciently strong, it may become possible for the quasiparti- 
cle states to delocalize, leading to a thermal metal phase. 
Consider, in particular, the situation where the ultimate 
destruction of the superconductivity at Hc2 leads to a 
normal metal. The normal metal is characterized by dif- 
fusive transport of charge, spin, and energy at the longest 
length scales. It is natural to expect that generically, the 
spin and thermal diffusion already exist in the supercon- 
ducting phase just below Hc2- In other words, we expect 
that a transition to a normal metal at Hc2 occurs from 
the superconducting thermal metal phase. Thus, in this 
case, there are three distinct zero temperature phases as 
the magnetic field is increased from just above Hd to 
just above Hc2- The superconducting thermal insulator 
at low fields (above Hd) first undergoes a delocalization 
transition to the superconducting thermal metal at some 
field Hc4 < Hc2 before the superconductor is destroyed 
to form the normal metal (see Figure 0). 
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FIG. 2. A possible phase diagram for a dirty Type // su- 
perconductor 



In the other case where the destruction of superconduc- 
tivity at Hc2 leads to a localized insulator, it is natural to 
expect that the transition occurs from a superconductor 
where the quasiparticles are localized, i.e directly from 
the superconducting thermal insulator (see Figure 0) . 
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FIG. 3. Another possible phase diagram for a dirty Type 
// superconductor 



III. SUPERCONDUCTING THERMAL METAL 

In this section, we examine the properties of the su- 
perconducting thermal metal phase in more detail. In 
this phase, there is diffusive transport of spin and en- 
ergy. The quasiparticle density of states at the Fermi 
energy is non-zero and finite. 

We may quantify the spin transport by defining a "spin 
conductivity" in analogy with the electrical conductivity 
for charge transport. The role of the chemical potential 
is played by a Zeeman magnetic field B coupling to, say, 
the z-component of the spin. The analog of the electric 



field is thus the spatial derivative of the Zeeman field. 
The spin conductivity Us thus measures the z component 
of the spin current j| induced in the system in response 
to an externally applied spatially varying Zeeman field 
along the z-direction of spin: 
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(13) 



Here g is the gyromagnetic ratio, and fis the Bohr mag- 
neton. It is easy to show that as satisfies an Einstein 
relation 



DsXo 



(14) 



where D^ is the spin diffusion constant, and xo is the spin 
susceptibility. In the approximation of ignoring quasipar- 
ticle interactions, xo is simply proportional to the quasi- 
particle density of states po a-t the Fermi energy: 



Xo = -. Po- 



(15) 



In the thermal metal phase, <Ts is finite and non-zero at 
zero temperature. 

In the vortex phase, the diffusion parallel to the mag- 
netic field, described by the diffusion constant Dy, is 
mainly along the core of the vortices, while the diffusion 
perpendicular to it, with diffusion constant D±, is due 
to in-plane motion between vortices. The latter depends 
strongly on the inter-vortex tunneling strength, and in 
general we expect the diffusion to be highly anisotropic. 
For ease of presentation, we will, for the time-being, as- 
sume that the diffusion is isotropic. When appropriate, 
we will take into account the anisotropic diffusion. 

The energy diffusion is measured by the more familiar 
thermal conductivity k. This too satisfies an Einstein 
relation: 



DtC, 



(16) 



where D^ is the thermal diffusion constant, and C the 
specific heat. In the approximation of ignoring quasipar- 
ticle interactions, the specific heat is determined by the 
density of states. In particular, in the limit T ^ 0, where 
T denotes the temperature, we have 



C = —k'sPoT. 



(17) 



Furthermore, in the non-interacting theory, as both spin 
and energy transport are by the quasiparticles, the cor- 
responding diffusion constants are the same: 



Ds = Di 



(18) 



This then leads to a "Weidemann- Franz" law relating the 
spin and thermal conductivities: 
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(19) 



Though the quasiparticle density of states is finite and 
non-zero at the Fermi energy, as we show below, it varies 
as \E on moving in energy [E) away from the Fermi en- 
ergy (See Fig. 0). This is in sharp contrast to a normal 
metal of non-interacting electrons , and has its origins 
in quantum interference effects specific to the supercon- 
ductor. This result can be established in a field theoretic 
analysis of quantum interference effects in the supercon- 
ducting thermal metal phase [pj. Here, we provide in- 
stead a simple physical explanation of the effect using 
a semiclassical argument developed earlier in a different 
context 101 . 




FIG. 4. Quasiparticle density of states of the superconduct- 
ing thermal metal 



The symmetry Eqn. of the Bogoliubov-deGennes 
Hamiltonian implies that the amplitude iGij^ap — 
{ia\e ~\jP)(){t) for a d-particle to go from point j, 
(pseudo)spin (3 , to point i, spin a, satisfies the relations: 



GijMit) =G*j-^||(i)- 
The Fourier transform of this amplitude is given by 

,. , 1 



(20) 
(21) 
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The quasiparticle density of states at an energy E away 
from the Fermi energy may be obtained from this in the 
usual manner: 

p[E) = -^/"i (G.^,TT {E + 11^) + (T^i)) , (22) 

where the overbar denotes an ensemble average over im- 
purity configurations. Consider now the return amplitude 
Gii^]]{t). This can be written as a sum over all possible 



paths for this event. Consider in particular the contri- 
bution from the special class of paths where the particle 
traverses some orbit and returns to the point i in time 
t/2 with spin down, and then traverses the same orbit 
again in the remaining time and returns with spin up. 



Using the symmetry relation Eqn,21, this contribution 
to iG'iiji(i) can be written as 



It is amusing to note that this correction has the same 
form as the Altshuler-Aronov effects in a diffusive, inter- 
acting normal metal, though the physics is quite different. 
Later on in the paper, when we consider interaction ef- 
fects , we will show the existence of an Altshuler-Aronov 
correction of the same form in the superconducting ther- 
mal metal as well. 
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Now |Gm^||(|)P is just the probability for the event 
i "f-^ J i in time t/2. For large t, this is half the to- 
tal return probability P{t). This in turn is determined 
by the condition that the d-particles diffuse through the 
system. For diffusing particles in three dimensions, the 
return probability is 
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(23) 



This then leads to an energy dependent correction to the 
density of states: 
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(24) 



If we take into account the anisotropic nature of the 
diffusion in the vortex phase, then the same result holds, 

but with an effective diffusion constant Dg = (D^Dy) ^ . 
The energy dependence of the density of states has im- 
portant consequences for the low temperature thermody- 
namics of the superconducting thermal metal phase. In 
particular, the specific heat at low temperature behaves 



as 
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where 7 — ^k^p^^ and the constant factor is given by 



6„,^iSl.^--i,)C(5/2). 
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Here C(5/2) — X^n —tj2- Note that hcq is a universal 
constant. Similarly, the spin susceptibility at low tem- 
peratures behaves as 
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The constant bs„ is again universal: 
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(27) 



(28) 



IV. SUPERCONDUCTING THERMAL 
INSULATOR 

We now consider the properties of the superconducting 
thermal insulator phase. By definition, this is a supercon- 
ductor where the thermal conductivity has the limiting 
form k/T ^ as T ^ 0. Similarly, the spin conductiv- 
ity is given by CTs = at zero temperature. Thus this 
phase is a superconductor for charge transport, but an 
insulator for thermal and spin transport. 

In contrast to conventional disordered insulators, the 
density of quasiparticle states actually vanishes in the su- 
perconducting thermal insulator. This can be seen by the 
following simple argument [H . Consider the Hamiltonian 
(^) in the limit of strong on-site randomness and weak 
hopping between sites. In the extreme limit of zero hop- 
ping, the sites are all decoupled. At each site, the Hamil- 
tonian in terms of the d-particles satisfies the SU{2) in- 
variance requirement ay Hay = —H*. This constrains 
the Hamiltonian to be of the form H — taz + Arax + J^^iay 
with t,Ar,Ai random. Physically, t can be thought of 
as a random on-site chemical potential, and Ar,Ai as 
the real and imaginary parts of the random on-site BCS 
order parameter A. Considering now the case where the 
joint probability distribution of t, A^, A.; has finite, non- 
zero weight at zero, we see immediately that the disorder 
averaged density of states vanishes as E^ . Now consider 
weak non-zero hopping. In the localized phase, perturba- 
tion theory in the hopping strength should converge, and 
we expect to recover the single site results at asymptot- 
ically low energies. If the joint probability distribution 
of t, Ar, Ai has vanishing weight at zero (as happens for 
instance, for A^, A^ non-random and non-zero, and only 
t random) , then the density of states only vanishes even 
faster. To get a density of states that vanishes slower 
than E'^ , or is finite at the Fermi energy requires a diverg- 
ing probability density at i = A^ = A^ = which is pre- 
sumably unphysical, and definitely non-generic. Thus, 
we conclude that the (disorder averaged) quasiparticle 
density of states vanishes in the superconducting ther- 
mal insulator phase in the absence of quasiparticle inter- 
actions. 

We now demonstrate the validity of this argument by 
direct numerical calculation of the density of eigenstates 
of the Bogoliubov-deGennes equations in the spin insu- 
lating phase. The simulations were done in one spatial 
dimension. The advantage of doing so is three-fold: as 



localization effects are strongest in one dimension, it is 
easier to access the properties of the localized phase in 
d = I. Further, it is possible to go to fairly large system 
sizes in d = 1 and hence the results are more reliable. Fi- 
nally, we expect that the properties of the localized phase 
are qualitatively the same in any dimension. Hence it is 
sufficient to consider the d = 1 case. A physical realiza- 
tion of a one dimensional system to which our results are 
directly applicable is obtained by considering the quasi- 
particle states in the core of a single, isolated vortex in 
the presence of disorder. 

To simulate the density of states of the d-particles, we 
employ the Hamiltonian in Eqn. |5|: 
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(29) 



where now the sites i and j reside on a one-dimensional 
lattice with periodic boundary conditions. It is conve- 
nient to picture the Hamiltonian in terms of coupled spin- 
up and spin-down sublattices. 

We now employ the Hamiltonian in Eq.E3 to explore 
the density of states numerically for various models and 
degrees of disorder. We begin with a model which shows 
a gap in the density of states in the absence of disorder. 
We set the nearest neighbour coupling tu^i to a constant, 
t, and take A = (5yA° to be on-site and real. The pure 
Hamiltonian can be diagonalized trivially. The resulting 
single particle excitations have a dispersion 



E = ±y/{A0)^+4:t^cosk, 



(30) 



with a gap 2 A" about E = 0, and a bandwidth 
2-\/(A°)2 + At"^ . We now introduce disorder by allow- 



(1) 



.(2) 



ing the on-site couplings tu, A^^ and A^^^"^ to take the 
values y^, A° -t- V'^ and V^ respectively. The V are ran- 
dom variables drawn from a uniform distribution with 
zero mean and width W which acts as a measure of the 
disorder strength. Note that a non-zero value of A^^) 
breaks time-reversal invariance, as one would expect for 
the vortex phase. 

As seen in FigJ|, impurities begin to fill in the gap. 
(The precise manner is specific to the distribution of dis- 
order, which we verified by using different forms for the 
probability distribution of the V"^). One can observe the 
symmetry about £^ = which is a result of the particle- 
hole symmetry of the Bogoliubov deGennes Hamiltonian. 
As the disorder strength is increased, there is an increas- 
ing density of states in the gap. But the density of states 
at the Fermi energy nevertheless always vanishes. Closer 
examination of the density of states near the Fermi en- 
ergy (see Figure H) shows that it actually vanishes as 



p{E)=A\E\^ 



(31) 



where A is a constant. This power law is exactly what is 
predicted by the simple argument outlined at the begin- 
ning of this section. 
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FIG. 5. Density of states in the superconducting thermal 
insulator showing evolution with increased disorder, W. The 
energy scale is set by the bandwidth of the pure dispersion in 
Eq.po[ 




FIG. 6. Density of states of the superconducting thermal 
insulator at low energies. Here, the constant 'A' of Eq.Bil de- 
creases with increasing disorder strength W that ranges from 
0.8 to 2. 



We now consider the situation where the density of 
states is a constant in the absence of disorder. We set 
the nearest neighbor coupling iii+i and gap A^^+i to be 
real constants t and A respectively. We then obtain for 
the dispersion of the single particle excitations the form 



£; = 2v/A2+t2cos/c, 



(32) 



with a bandwidth 4\/A2 + t'^. We introduce disorder by 
letting the on-site coupUngs tu , A^-^ ' and A^-^ ' take the 
values V^, V^ and V'^ respectively, where the V^'s are 
random variables as specified for the case with the pure 
gapped dispersion. As seen in Fig. m, disorder reduces 
the density of states at the Fermi energy, ultimately forc- 
ing it to vanish as E^. 
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FIG. 7. Evolution of the density of states with disorder ,W, 
for a model where the density of states is constant in the clean 
limit. Here the energy scale is set by the bandwidth specified 
by the pure dispersion of EqJ32. 



In summary, the superconducting thermal insulator 
phase has the remarkable feature that the quasiparticle 
density of states actually vanishes at the Fermi energy. 
This is in striking contrast to a conventional Anderson 
insulator, and has several obvious consequences for the 
low temperature thermodynamic properties of the phase. 

The vanishing density of states also has consequences 
for thermal and spin transport at finite temperature 
which is presumably through variable-range hopping. 
This can be seen as follows: for a hop between two local- 
ized quasiparticle states separated by a distance i?, the 

_R 

overlap of the two states ^ e « where ^ is the localiza- 
tion length. The typical energy separation En between 
these states is determined by the density of states. When 
this vanishes as E^, the total number of states in a radius 
R in an energy interval E about the Fermi energy is of 
order E^B? (in three dimensions). Therefore the typical 



spacing Er 



The total rate for hopping a distance R 



is proportional toe '■ e ~^ rt) _ xhe first term in the expo- 
nential is the square of the wave- function overlap, and the 
second is the activation energy cost for the hop. (Here c 
is some unknown constant determined by the density of 



states). The total hopping rate, as measured by the spin 
or thermal conductivities, is obtained by summing over 
hops of all possible distances R. For low temperature, it 
is clear that the sum will be dominated by hops of size 
R ^ -j=. Thus the total hopping rate (and hence the 

spin/thermal conductivities) '-^ e~ V It^J. Note that the 
exponent 1/2 differs from the conventional Mott expo- 
nent of 1/4, and is due to the density of states vanishing 
asE"^. 



V. INTERACTIONS 

Our discussion of the properties of the two super- 
conducting phases has thus far been based on the non- 
interacting quasiparticle Hamiltonian of Eqn. |l| or Eqn. 
y. In this section, we consider the effects of interactions 
between the quasiparticles. A discussion of the micro- 
scopic origin of these interactions is provided in Appendix 
|A|. Here we take a more phenomenological approach. It is 
important to note that the interactions are short-ranged 
in space - the long-range Coulomb repulsion between the 
underlying electrons is screened out by the condensate. 
We keep only the interactions in the triplet channel to 
get the Hamiltonian 



H ^Tio + Hint , 

T-iint = I d'^xutix - y)S{x) ■ S{y), 



(33) 
(34) 



where Ut is short-ranged, and S{x) is the spin density. 
Note that as the charge density is not a hydrodynamic 
mode, the singlet interaction is expected to be quite in- 
nocuous, and we simply drop it. As discussed in Ap- 
pendix H, the sign of ut could be either positive or neg- 
ative, with negative ut corresponding to repulsive inter- 
actions. 

Interactions can also be included instead in the lattice 
model Eqn. y. We expect that the universal properties 
of the two superconducting phases are insensitive to the 
detailed form of the interaction Hamiltonian. So we may 
consider any short-ranged interaction with the right sym- 
metries (spin conservation) . A particularly simple choice 
is provided by an on-site Hubbard interaction: 






1), 



(35) 
(36) 



where rii = c^Ci is the number operator for site 'i'. For 
repulsive interactions, U > 0. 

In the rest of this section, we will use these model in- 
teraction Hamiltonians to discuss their effects on the two 
superconducting phases. 



A. Superconducting thermal metal 



Clearly then, we have 



Interaction effects lead to significant changes in the 
properties of diffusive normal metals Q , as was shown by 
Altshuler and Aronov. In the superconducting thermal 
metal phase, it is natural to expect similar effects due 
to the interplay of spin diffusion and interactions. We 
demonstrate this below with some simple calculations. 

Consider the free energy of the interacting quasiparti- 
cles in the superconducting thermal metal. To calculate 
this, it is convenient to pass to a functional integral rep- 
resentation for the partition function: 



Z= \ [Vijj]e 



-(So+S„t) 



(37) 



d 



So^ J drd^x (V(a;, r)— V(a;, r) + 'Hoii', ^)), (38) 

S^nt - f drd'^xd'^y mix - y)S{x, r) • Siy, r) (39) 

As we show below, all the singular corrections to the free 
energy come from the diffusive nature of the spin fluctu- 
ations. (The charge density does not diffuse, and hence 
the singlet interaction does not contribute to any singu- 
lar corrections ;we are justified in dropping it). The first 
order correction to the free energy due to interactions is 
given by 

VpAT = (S.nt) 

drd'^xd'^y Ut{x~y){S{x,T)-Siy,T)), (40) 



where V is the system volume, /3 the inverse tempera- 
ture, ^ the free energy density, and AT the correction 
to the free energy density due to interactions. The ex- 
pectation value on the right hand side is to be evaluated 
in the non-interacting theory. The overline denotes an 
average over all realizations of the disorder. The expec- 
tation value can be related directly to the spectral weight 
for spin density fluctuations as follows: 

{S{x,T)-S{y,T))^{S{x,0)-S{y,0)) 



E 



= -/3S„ 



(n|S'(a;)|m) • (■m\S{y)\n) 



where |n),|m) are exact eigenstates (with energies 
En , Em respectively) of the non-interacting Hamiltonian 
Ho, and Zq = Tre^^^" is the corresponding partition 
function. Now the spectral weight for spin density fluc- 
tuations in the non-interacting theory can be expressed 
as 



X"{x,y;uj) =^ 



-/3B„ 



Zo 



(n|S'(a;)|m) • {■m\S{y)\n) 



(2i)J(„ - '^- ^■> ) (1 - e-'"") 



(5(.,0).s„,o)).£^f:<:^^ (41) 

Upon averaging over the disorder, translation invariance 
is restored, and we have 



X"(x,y;uj) = 



d'^q 



(2^) 



-eu^-y)x {q,oj). 



(42) 



The singular contributions to the free energy all come 
from the small q, lo behaviour of x ■ This is entirely de- 
termined by the diffusive nature of spin transport. 



X (9,^) 



GhasUjq^ 



OJ^ 



(43) 



where as is the spin conductivity and Ds is the spin dif- 
fusion constant. (There is an extra factor of 3 in the ex- 
pression above from the standard form , which is due to 
summing over all three components of the spin) . We can 
now evaluate the free energy correction (See Appendix H 
for details). The result is given by 



AT{T) ~ AJ^(O) utT^/^ 



(44) 



in three dimensions. (Here we have ut = J d'^xut{x)). 
This then leads to a correction to the specific heat due 
to interactions of the form 



AC = br 



\2TThDs 



3/2 



where the coefficient bd is given by 



bci = --^ut 



Xo 



8 {9^^B)' 



rfcsC(5/2). 



(45) 



(46) 



Note that this is similar in form to the correction due to 



quantum interference pointed out in Section [II, but has a 
different physical origin. (As before, with anisotropic dif- 
fusion, Dg is replaced by the effective diffusion constant 

(yD\D\\) ^). The total specific heat in the superconduct- 
ing thermal metal phase therefore is 



C^-fT + b, 



ksT 
2TrhDs 



3/2 



(47) 



where the constant be has contributions from both quan- 
tum interference and interaction effects: 



^cq 



(48) 



Note that b^q is a universal number, while b, 



depends on 

the interaction strength and the zero temperature spin 
susceptibility. This should allow experimental determi- 
nation of the importance of interaction effects by simul- 
taneous measurements of the specific heat and the spin 



or thermal diffusion constant. The deviation of the coef- 
ficient be from the universal value bcq gives a measure of 
the interaction strength. 

Similar effects exist in the spin susceptibility as well. 
To see this, consider evaluating the free energy in the 
presence of an externally applied Zeeman magnetic field 
coupling to the z component of the spin. This does not 
affect the diffusion of Sz , but the diffusion of Sx , Sy gets 
cut off. The precise manner in which this happens can 
be found from hydrodynamic considerations by including 
the effect of the precession of the spin density under the 
external magnetic field in the classical diffusion equa- 
tion. To that end, first note that in the presence of a 
small magnetic field B coupling to the z-component of 
the spin, the ground state has a spin density 



S 



- Xq „ 
z B. 

gtJ-B 



(49) 



Consider small deviations of the local spin density from 
the ground state: 



S{x, t) = ( 5Sx{x, t),5Sy{x, t), -^B + 5S,Xx, t) 

9f^B 



(50) 



Quite generally, this satisfies the equation of motion 



as 

dt 



+ Sx [gtJLBB) = -ddl, 



(51) 



where i = 1, 2, 3 is a spatial coordinate index and jl is 
the spin current vector in the spatial direction i. The 
second term on the left hand side arises from the pre- 
cession of the spin under the external magnetic field. In 
the absence of this term, the equation above reduces to 
the familiar continuity equation expressing spin conser- 
vation. To derive the 5*^5" correlator, introduce a small 
additional magnetic field b{x,t) coupling to Sx- Then, 
the spin currents j",i|'* are related to gradients of the 
spin density and the field b through 



rs 



= ~D.,d,.ssy, 



(52) 
(53) 



where Ds , as are the spin diffusion constant and spin 
conductivity respectively. We may now determine the 
response of the spin density to the field b from these 
equations: 



<5s+(fc,w) 



Ds 



D,k'^ -iB 



-i{uj + B) + Dsk'^ 



b{k,u), (54) 



where B — gusB. We have used the Einstein relation 
EqnJIj to express xo '^^ terms of (Tg and Dg. A similar 
expression, but with B — > — i?, holds for Ss~{k,u!). The 
response function Xxx{k, uj) can now be read off to be 



Xxxik,uj) 



2Ds 



-i[uj + B) + Dsk'^ 



(B ^ -B) 



(55) 



This then determines Xxxi^T^) through the equality 



Xxx{k,^) = 2hlmxxx{k,uj). 
The result is 

Xxxiq,'^) = TT-^si^k"^ 



(56) 



1 



{co + By + {Dsk^ 



(B 



-B) 

(57) 

Exactly the same expression holds for x . Notice that 
the form of the spectral weight given in Eq.H3 may be 
obtained from the above by setting the magnetic field B 
to zero. 

It is now possible to evaluate the change in free en- 
ergy due to the magnetic field. For details, see Appendix 
y|. The susceptibility is then obtained by differentiating 
with respect to the field. We find, for the interaction 
correction, 



Axo(T) 



keT 
D, \2TThD. 



1/2 



(58) 



in three spatial dimensions. The coefficient bsi is given 

by 



8^3/2 ' 



(59) 



where T = f°° ^A. ( ^i—\ 

The full low temperature behaviour of the spin suscep- 
tibility is then given by 



X{T) - Xo = 



keT 
D, V 2'KhD, 



1/2 



(60) 



where again bg has contributions from both quantum in- 
terference and interaction corrections: 



(61) 



Note again that bgq is a universal constant while bsi de- 
pends on the interaction strength. (With anisotropic dif- 

fusion, Ds is replaced by (D'^D||) ^ ). Thus measurements 
of the temperature dependence of the spin susceptibility 
can also be used to infer the relative importance of inter- 
action and quantum interference effects. 

The simple calculations above demonstrate the exis- 
tence of Altshuler-Aronov effects in the thermodynamic 
properties of the superconducting thermal metal pO[ . It 
is also possible to show that interaction effects lead to 
a suppression of the electron tunneling density of states. 
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For the tunneling of a d particle, the result can be es- 
tablished in quite a straighforward way along the lines 
of Ref. ||T^. However, the physical tunneling process in- 
volves tunneling of an electron. In Appendix pi we show 
that under certain further approximations, the tunnel- 
ing density of states of electrons is the same as that of 
the d particles. Thus, interaction effects lead to a sup- 
pression of the physical tunneling density of states in the 
superconducting thermal metal as well. 



B. Superconducting thermal insulator 

We now move on to consider the effects of interactions 
in the superconducting thermal insulator. In a phase 
with a gap in the quasiparticle spectrum, weak inter- 
actions are irrelevant and lead to no significant effects. 
(This is analogous to the insignificance of weak interac- 
tions in an ordinary band insulator). So we consider the 
more interesting case of a gapless superconducting ther- 
mal insulator . We argued in Section |^ that, in the 
absence of interactions, the disorder averaged density of 
states vanishes, generically as |i?P, on approaching the 
Fermi energy. Despite this, we show here that arbitarily 
weak repulsive interactions lead to the formation of free 
paramagnetic moments. This result is quite analogous 
to what happens in a conventional Anderson insulator; 
arbitrarily weak repulsive short-ranged interactions lead 
to the formation of free paramagnetic moments. (Note 
however that the density of states is a constant in the 
conventional Anderson insulator) 

The discussion is simplest in terms of the lattice 
Hubbard-type Hamiltonian Eqn. p9. As in Section pv| , 
consider the case of strong on-site randomness and weak 
hopping. In the limit of zero hopping, the resulting single 
site problem can be solved exactly even in the presence 
of the interaction term. To that end, consider the single 
site Hamiltonian, 









a44 



CiCT, 



Hu = t/n(n-l). 



(62) 
(63) 

(64) 



with n = ^^ c)^Ca- For the time being, we assume [/ > 0, 
making the interactions repulsive. It is convenient to go 
to the d representation: 



"Ho = d^ [ta^ + ArCTx + ^iO-y) d, 
Hu = U[{d^a,df+d^cr,d] . 



(65) 
(66) 



Here A^^i are the real and imaginary parts of A. We 
will assume that t,Ar,Ai are random variables with a 
probability distribution that has finite non-zero weight 
when all three variables are zero. As we argued in Sec- 
tion III, in that case, the quasiparticle density of states 
of the non-interacting Hamiltonian goes to zero as \E\'^. 



To diagonalize the full interacting Hamilonian, note 
that the d particle number is conserved even with inter- 
actions. Thus we may look for eigenstates with fixed d 
particle number n^ = X^a'^a'^"- T^^ physical spin is 
determined entirely by Ud (Eqn. ^ through the rela- 
tion 5*2 = ^{ud — !)■ For the single site problem, the 
Hilbert space consists of four states: |0), | t) = '^t|0), 

I i) = 4l0)' ^^d I ti> = 44 10)- The states |0) and | Ti) 
are immediately seen to be eigenstates with energy 0. Di- 
agonalizing the Hamiltonian in the space of the remaining 
two states | t), | i), we find two other eigenstates |-|-), |— ) 
with the eigenvalues E± = U ± y^it + Uy + jAp. Note 
that the lower of the two eigenvalues £'_ is negative if 
—2tU < t^ + |Ap. (The other eigenvalue Ej^ is always 
positive). If -E- < 0, then the ground state is the state 
|— ). This lies in the subspace with Ud = 1^ and hence has 
physical spin — 0. If E^ > 0, then the ground state is 
two fold degenerate with both |0), and | t|) having zero 
energy. These correspond to states where the ground 
state has a physical spin 1/2. 

Thus a free magnetic moment is formed in the ground 
state of a single site whenever the condition t^ + |Ap < 
~2tU is satisfied. For a collection of sites with a generic 
random distribution P{t,ArAi) of i, A^, and A^, there 
will always be some weight where this condition is satis- 
fied. Hence there will be a finite probability of forming a 
free moment at any site. This leads to a finite density of 
magnetic moments for the full system. 

It is natural to expect that inclusion of weak hopping 
between sites does not change the result above, so long as 
we are in the localized phase. So we conclude that arbi- 
trary weak repulsive interactions lead to the formation of 
free paramagnetic moments in the superconducting ther- 
mal insulator . 

What is the ultimate fate of these magnetic moments 
at low energies? There are two generic possibilities: the 
spins can freeze into a spin glass phase, or stay unfrozen 
in a phase with random singlet bonds |12[ between pairs 
of spins. It is a difficult matter to decide on the condi- 
tions for realizing either of these possibilities, and we will 
not attempt it here. 

The analysis above has assumed that the interaction 
was repulsive. If the interaction is attractive, i.e U < 0, 
the solution of the single site problem shows that the 
ground state has no net spin. Thus there is no local 
moment instability in that case. The spin susceptibility 
vanishes at zero temperature, as in the non-interacting 
case. 



VI. SYSTEMS WITH TIME REVERSAL 
INVARIANCE 



In this section, we will briefly consider systems with 
time reversal invariance. An s-wave superconductor with 



11 



weak or moderate impurity scattering has a gap in the 
quasiparticle excitation spectrum, and hence is in the su- 
perconducting thermal insulator phase. Strong impurity 
scattering will destroy the superconductor. When this 
happens, the resulting phase is most likely an insulator. 
However we see no reason of principle forbidding a transi- 
tion to a normal metal (in three spatial dimensions) when 
the superconductor is destroyed. We propose that the 
transition then occurs from the superconducting thermal 
metal phase, i.e as the impurity strength is increased, 
there is a transition from the superconducting thermal 
insulator to the superconducting thermal metal which 
precedes the ultimate transition to the normal metal (See 
Figure g). On the other hand, if the transition is to the 
insulator, we propose that it is directly from the super- 
conducting thermal insulator (Figure y). 



SUPERCONDUCTING 
THERMAL EVSULATOR 



SUPERCONDUCTING 
THERMAL METAL 



NORMAL METAL 



W 



FIG. 8. Zero temperature phase diagram for the supercon- 
ductor-normal metal transition in the presence of time rever- 
sal invariance. The parameter VF is a measure of the strength 
of the disorder. 
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FIG. 10. Density of states in the T -invariant supercon- 
ducting thermal insulator showing evolution with increased 
disorder, W. The energy scale is set by the bandwidth of the 
pure dispersion of Eqn. ku of Section 



The density of states clearly vanishes at zero energy. 
Closer examination of the low energy behaviour (See Fig- 
ure |ll|) shows that it actually vanishes as \E\ in agree- 
ment with the simple argument above. 



FIG. 9. Zero temperature phase diagram for the supercon- 
ductor-insulator transition in the presence of time reversal 



invariance. 



The properties of the superconducting thermal metal 
phase in this case (with time reversal symmetry) are quite 
similar to the case where time reversal symmetry is ab- 
sent. Differences exist, however, in the superconducting 
thermal insulator phase. The heuristic argument 0], in 
the beginning of Section [V; for the quasiparticle den- 
sity of states now shows that it vanishes in this case as 
well, but only as fast as or faster than |_E|. As in Sec- 
tion p^, we will provide supporting numerical evidence 
for this statement by calculating the density of states of 
the Bogoliuobov-deGennnes equations appropriate for a 
time reversal invariant system in one dimension. In or- 
der to preserve time reversal symmetry, we set the imag- 
inary part of the gap function A*-^) in the Hamiltonian of 
Eqn|23 of Section f^ to zero. Apart from this important 
difference, the models that we use here are the same as 
those of Section |^. Again, we first display results for a 
one dimensional model where there is a gap in the clean 
limit. 




Energy 



FIG. 11. Density of states in the time reversal -invariant 
superconducting thermal insulator at low energies. Decreas- 
ing slope corresponds to increasing disorder strength, W, that 
ranges from 0.8 to 4. 



We also considered the situation where the density of 
states is a constant in the clean limit. The results are 
shown in Figure O. Again, the density of states vanishes 
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linearly for strong disorder. 



C 



\H-Hr. 



(67) 



1 . Pure system 





Energy 





FIG. 12. Evolution of density of states with disorder, W, 
in model where it is constant in the clean limit. Here, the 
energy scale is set by the bandwidth of the pure dispersion 
given by Eqn.|| of Section ^ 



VII. PHASE TRANSITIONS 



In this section, we study some general properties of 
the various transitions between the phases that we have 
discussed (see Fig.|l^. The most novel phase transition 
in our phase diagram is the one between the supercon- 
ducting thermal metal and the superconducting thermal 
insulator. This is a localization transition that occurs 
inside the superconducting phase, and should be accessi- 
ble experimentally. On approaching the transition from 
the thermal metal side, the spin conductivity goes to 
zero continuously. Similarly, the low temperature ra- 
tio A = k/T, which is non-zero in the thermal metal, 
vanishes at the transition. The low temperature ther- 
mal transport in the superconducting thermal insula- 
tor is presumably through variable-range hopping, and 
K r^ e ^ "^ ' with X = 1/2. In the critical region near 
the transition point, various physical quantities are ex- 
pected to have universal singular behaviour. In the model 
of non- interacting quasiparticles, it is possible to develop 
systematic calculations of the universal critical exponents 
p3[. However, we expect universal properties even with 
interactions present. 

On approaching the transition at zero temperature by 
varying the field strength towards Hd, there is a length 
scale ^ (which may be interpreted as the quasiparticle 
localization length in the localized phase) that diverges 
as 



Similarly, moving away from the critical point by turn- 
ing on a finite temperature introduces a length scale ^t 
which behaves as 
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T" 



(68) 



Note that the two equations (p^) and (BSl) define the two 
critical exponents v and z. 

On approaching the transition point from the delocal- 
ized side, the coefficient A introduced above goes to zero. 
Precisely at the transition point H = -ffc4, but at finite 
temperature, k{T) ^ 7^1+0 with > being a univer- 
sal exponent. In general, for fields close to Hd and low 
temperatures, we may write down a scaling form 



^ = c^T^Y (I 



(69) 



The constant ci is non-universal, while the function Y is 
universal. Equivalcntly, we may use Eqns. ( |67| ) and ( |68| ) 
to write 



K 

T 



Cl 



j^4>Y 



\H~H. 



C2- 



c4| 



T± 



(70) 



Here C2 is also non-universal, and F is a universal scal- 
ing function such that 1^(0) is a finite constant. Further, 
requiring that k ~ T as T — )■ on the metallic side, we 
see that Y{x -^ +00) ~ x^'^'^. This then implies that the 
coefficient A introduced above vanishes, on approaching 
Hc4, as (H — HcaY'^'^ ■ Further, if we make the assump- 
tion that the fixed point theory controlling the transition 
obeys hyperscaling, then conventional scaling arguments 
can be used to show the exponent equality 



d-2 



(71) 



Similar scaling forms, but with different exponents and 
scaling functions, describe the transition in the time re- 
versal invariant case as well. 

Our phase diagram has important implications for the 
phase transitions at zero temperature where the super- 
conducting phase is destroyed. For the superconductor- 
normal metal transition, we suggest that the transition 
is generically from the superconducting thermal metal 
phase, i.e, from a superconductor with spin and en- 
ergy diffusion at T = 0. Further this implies that 
the superconductor-normal transition generically occurs 
from a gapless superconductor with a finite, non-zero 
density of quasiparticle states at the Fermi energy. The 
presence of gapless quasiparticle excitations in the super- 
conducting side should affect strongly the critical proper- 
ties of the SC-normal metal transition: previous theoret- 
ical treatments of this transition p4[ | which have ignored 
this feature are hence expected to be incorrect. Simi- 
larly, the superconductor-insulator transition is generi- 
cally from the superconducting thermal insulator. 



13 



NORMAL METAL 



SUPERCONDUCTING 
THERMAL METAL 



NORMAL INSULATOR 



SUPERCONDUCTING 
THERMAL INSULATOR 



FIG. 13. Zero temperature phase diagram showing transi- 
tions where the superconducting phase is destroyed 



As shown in Fig. n3, direct transitions from the su- 
perconducting thermal metal to the normal insulator or 
from the superconducting thermal insulator to the nor- 
mal metal are possible at a special multicritical point. 
We will however not attempt to describe the properties 
of such a multicritical point here. 



VIII. DISCUSSION 

One of the main purposes of this paper is to point out 
that all superconductors fall into one of two categories 
- those that, apart from being superconducting, share 
many properties with conventional metals, and those that 
share many properties with conventional insulators. This 
distinction can be made quite precise, and indeed corre- 
sponds to the existence of two distinct zero temperature 
superconducting phases which are distinguished by the 
nature of quasiparticle transport. In this final section, 
we will discuss various experimental implications. This 
will be followed by a discussion of various peripheral mat- 
ters that have been omitted from our considerations so 
far, and their effects on experimental systems. 



A. Experiments 

A powerful way of probing quasiparticle transport in 
a superconductor is through thermal conductivity mea- 
surements. In the superconducting "insulator" phase, 
the ratio n/T goes to zero as the temperature goes to 
zero. On the other hand, in the superconducting "metal" , 
k/T goes to a constant as the temperature goes to zero. 

The Type // s -wave superconductor offers a definite 
opportunity for tuning between these two phases. At low 
fields, this is in the superconducting "insulator" phase. 
Consequently, k/T ^ as T ^ 0. Upon increasing the 
field, there could, under conditions we have outlined, be a 
delocalization transition to the superconducting "metal" 
phase with k/T going to a constant. It should also be 
possible to explore the properties of the phase transition 
between the superconducting "insulator" and the super- 
conducting "metal" . Right at the critical point, the ther- 
mal conductivity k ~ T^+'l> with (/> > 0. We are not aware 



of any experimental investigations of this novel "metal- 
insulator" transition inside the superconducting phase so 
far. In performing the heat transport measurements, it 
is, of course, necessary to ensure that the phonon contri- 
butions have been subtracted out ||l5| . 

Another interesting experimental possibility is pro- 
vided by three dimensional dirty s-wave superconductors 
in the absence of any external magnetic fields. Upon 
varying the impurity concentration to destroy supercon- 
ductivity, if the transition is to a normal metal, we 
have proposed that it occurs from the superconducting 
"metal" phase. As the superconductor is in the super- 
conducting "insulator" phase at low impurity concentra- 
tions, there will then be a phase transition to the super- 
conducting "metal" with increasing impurity concentra- 
tion before the destruction of superconductivity. This 
too can be probed by thermal transport measurements. 
The phase transition in this case is similar to the "metal- 
insulator" transition discussed above, but belongs to a 
different universality class due to the presence of time 
reversal symmetry. 

Recent measurements [M of the low temperature in- 
plane thermal conductivity in the high temperature su- 
perconductors show that k/T goes to a constant as 
T ^ 0. Whether this is indicative of a true three dimen- 
sional superconducting "metal" phase stabilized by in- 
terlayer quasiparticle hopping, or just a two dimensional 
weakly localized superconducting "insulator" (see below) 
is difficult to ascertain at present. 

A number of experimental predictions follow from our 
study of the properties of either phase. We start with 
the superconducting "metal" . For systems with negligi- 
ble spin-orbit scattering, this phase is characterized by 
the presence of spin diffusion at zero temperature. Mea- 
surements of the spin diffusion constant should then show 
a non-zero value at zero temperature. The spin suscep- 
tibility, as measured, for instance, by the Knight shift is 
also predicted to saturate to a finite, non-zero value at 
zero temperature. Note however that due to quantum 
interference, and interaction effects, we predict a vT 
dependence in the temperature dependence of the spin 
susceptibility at low temperatures. Similarly the specific 
heat has a T^' ^ correction. As we emphasized in Section 
M, either the specific heat or the susceptibility measure- 
ments may be used to quantify the relative importance 
of interaction and quantum interference effects. 

In the superconducting thermal insulator phase, the 
spin conductivity is zero at zero temperature. Further, 
if the interactions are weak, then the vanishing density 
of states in the non-interacting quasiparticle model mani- 
fests itself in all thermodynamic properties. For instance, 
the specific heat vanishes faster than T"+^ with a — 1,2 
depending on whether time reversal is a good symme- 
try or not. Similarly, the spin susceptibility vanishes 
faster than T". Interaction effects become important at 
the lowest temperatures - formation of free paramagnetic 
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moments would initially give rise to a Curie term in the 
susceptibility which would then be altered at even lower 
temperatures due to exchange between these moments. 

So far, we have neglected the Zeeman coupling and 
spin orbit interactions. In making contact with experi- 
ments, it is essential to have some understanding of the 
effects of including these. We discuss that next. 



value of k/T as T ^ 0. In contrast to systems with con- 
served spin, the density of states in the non-interacting 
quasiparticle theory has a y/E enhancement at low 
energies. Similarly, in the localized insulator, in the non- 
interacting theory, spin-orbit effects allow Q for the pos- 
sibility of a finite density of quasiparticle states at the 
Fermi energy. 



B. Zeeman coupling 

Weak Zeeman coupling does not affect the existence 
of the two kinds of superconducting phases, but modi- 
fies their properties. For a magnetic field along, the z- 
direction of spin, the components of spin along the x, y 
axes are no longer conserved. This cuts off their dif- 
fusion at long length and time scales in the supercon- 
ducting "metal" phase. The z-component of the spin 
and the energy continue to diffuse though. Further, both 
the quantum interference and Altshuler-Aronov interac- 
tion corrections to the thermodynamic quantities in the 
"metallic" phase are cut-off by the finite Zeeman cou- 
pling. For instance, this leads to a field dependent spin 
susceptibility at the lowest temperatures of the form 



x{B)-xo-B^/\ 



(72) 



For the interaction correction, this is demonstrated in 
the Appendix. For the quantum interference contribu- 
tion, this result may be understood by noting that the 
Zeeman magnetic field acts as a "chemical potential" for 
the d-particles (as it couples to the physical spin ~ d- 
particle number). The density of states at the Fermi 
energy (and hence the spin susceptibility) is then given 
by Eqn. gj to be ~ B'^/'^ . On the spin-insulating side, 
in the approximation of ignoring interactions, the effect 
of finite Zeeman coupling depends on whether or not a 
hard gap exists in the quasiparticle excitation spectrum. 
If gapped, weak Zeeman coupling is innocuous, and can 
be ignored. On the other hand, if the system is gapless 
with a density of states vanishing as \E\°' (a — 1,2 de- 
pending on whether time-reversal is a good symmetry or 
not), then weak Zeeman coupling leads to a finite density 
of states at the Fermi energy, proportional to \B\°' . 



C. Spin-orbit effects 

Spin-orbit scattering, like Zeeman coupling, does not 
affect the existence of the two kinds of superconducting 
phases in three dimensions, but modifies their properties. 
In the presence of spin-orbit scattering, no component of 
the spin diffuses in the superconductor with delocalized 
quasiparticles. Energy continues to diffuse in this phase 
however. Thus this phase is characterized by a finite 



D. Two dimensional systems 

Through out this paper, we have focused on three di- 
mensional systems. Here we make some brief remarks 
about two dimensional systems. For systems with con- 
served spin, quantum interference effects lead [Ej to an 
absence of spin and energy diffusion at the longest length 
scales. The superconducting thermal metal therefore 
does not exist in two dimensions, at least if quasiparticle 
interactions arc ignored. The superconductor-insulator 
transition in two dimensions therefore occurs from the 
superconducting thermal insulator phase. Spin-orbit 
scattering can stablilize [plj^ a phase with delocalized 
quasiparticles in the two dimensional superconductor in 
the approximation of ignoring quasiparticle interactions. 
The resulting phase has a divergent k/T as T — > 0, and 
a divergent density of quasiparticle states at the Fermi 
energy 0. 

This research was supported by NSF Grants DMR-97- 
04005, DMR95-28578 and PHY94-07194. 



APPENDIX A: INTERACTION HAMILTONIAN 

Here, we provide some microscopic justification for the 
interaction Hamiltonian Eqn. B4. For simplicity, we con- 
sider a clean system in d = 2. This may be described by 
the model Hamiltonian 



H — Ho + Hcoul + Hel-ph, 

Ticoui ^Y. I ^l{x)Mx)y{x - x')iPlix')iP^,{x'). (A3) 



(Al) 
(A2) 



Here V{x — x') ~ \x-x'\ ^^ ^^'^ Coulomb interaction. 
Tiei-ph is the electron-phonon interaction which we do 
not specify in detail other than to assume that it leads 
to an effective attractive interaction at energies smaller 
than the Debye frequency ujd << Ep about the Fermi 
energy. We first imagine integrating out all modes p^ ] 
except those within ujo about the Fermi energy. In 
the resulting low energy theory, only a small number 
of qualitatively different interaction amplitudes are al- 
lowed due to geometrical restrictions imposed on the 
scattering processes p^. These correspond to a charge 
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density-charge density ("singlet" channel), spin density- 
spin density ("triplet" channel) interactions and to two 
interactions between spin singlet and spin triplet Cooper 
pair operators. Perturbative renornialization group argu- 
ments |l^ show that the singlet and triplet amplitudes 
are marginal up to one loop, while the Cooper channel 
amplitudes are marginally relevant if attractive. We as- 
sume such an attractive interaction leads to a flow to- 
wards a spin-singlet BCS superconductor. Treating the 
attractive spin-singlet Cooper interaction in mean field 
theory, and ignoring the other interactions in the sin- 
glet and triplet channels is equivalent to conventional 
BCS mean field theory. Going beyond the BCS theory 
requires reinstating the singlet and triplet interactions, 
and including fluctuations of the BCS order parameter. 
The order parameter fluctuations may be integrated out 
in the superconducting phase leading to an effective four 
fermion interaction which renormalizes the singlet ampli- 
tude. This then leads to an effective action for the quasi- 
particles in the superconductor which includes quasipar- 
ticle interactions. 

Consider now the system we have focused on the most 
- the Type // superconductor in a field in the presence 
of disorder. In principle, the interaction Hamiltonian 
contains both the singlet (charge-density) interactions 
and the triplet (spin-density) interactions. However, the 
charge density is not a hydrodynamic mode in the su- 
perconductor, and does not diffuse. Consequently, the 
singlet interaction is expected to be quite innocuous. We 
therefore retain only the spin triplet interaction. 

The sign of the (bare) triplet interaction ut is deter- 
mined by the balance between the repulsive Coulomb in- 
teraction amplitude in the triplet channel and the attrac- 
tive (phonon) interaction amplitude in the same chan- 
nel. This is, in principle, different from the balance in 
the Cooper channel where a net attractive interaction is 
required for superconductivity. Thus, Ut can be either 
positive or negative. It can be shown that ut < corre- 
sponds to repulsive interactions. 



APPENDIX B: ALTSHULER-ARONOV 
CORRECTIONS 

In this appendix, we provide details of the calculation 
of the singular corrections to the properties of the super- 
conducting thermal metal phase due to interactions. 



where Ut{k) = J d^x e~'^^''^Ut{x) is the Fourier trans- 
form of the triplet interaction Ut{x). The k integral is 
restricted to |fc| < A ^ j- where l^ is the elastic mean 

free path. (We have denoted dk = -^^)- The oj inte- 
gral runs from — oo to -|-oo. If the range of the short- 
ranged interaction ut{x) is much smaller than the mean 
free path (which we assume) , then the Fourier transform 
ut{k) may be approximated by it's value at /c = 0, i.e 
ut{k) K, Ut = J d'^xut{x). Subtracting out the zero tem- 
perature correction to the free energy, we get 



StAT = AT{T) - AT{0), 
5tAJ^ 



QTiutUs 



dk- 



^duj / 1 
2tt \ 1 - e-l^'^'^ 

/■oo 

dkk / duj 
Jo 



e{oj) 



LOp 



W2 + i:)2fc4 



Dlk^' 



In going to the second line, we have integrated over 
the angular coordinates of the momentum fc, and have 
reduced the frequency integral to one over positive uj 
alone. We now make the change of variables y = PTllo, 
X = y/(3hDsk to get 



tt^6t{AT) _ fkeT 
3hutO's \fiDs 



x,y 



V 



ev 



1 



Here /^ ^ = /p"' dx ^^ dy and xo = l^sj^- The inte- 
gral over X, y above can clearly be rewritten as 



6 



dx 



dy 



ev -1 



1 



y 



The first term contributes 0{T^) to the free energy, and 
is hence non-singular. All the singular corrections come 
from the second term. As the x-integral is ultraviolet 
convergent in this term, we set xq = oo and evaluate 
it explicitly to obtain the singular correction to the free 
energy: 



AT{T) - Ajc-(O) 



'dfiUtUs 

'4(27r)3/' 



-^^'^^m-T 



(Bl) 



Expressing CTs in terms of Dg and xo using the Einstein 
relation Eqn. |4[ and differentiating with respect to T to 
get the specific heat, we get the result quoted in Section 

0- 



1. Specific heat correction 



2. Susceptibility correction 



As shown in Section the correction to the free en- 
ergy at zero magnetic field is given by 



AT = &hc7, fdk^utik) (- 

J 27r \^-e 



-f3huj cj^ + D^k'^' 



The field dependent term in the correction to the free 
energy is 



AJ^{B) = 2hiit / dk 



^dLJ 1 



:Xxxik,uj). 



16 



The factor of two in front accounts for the contribution 
from both the SxSx and SySy correlators, and Xxxi^j'^) 
is given by Eqn. [5^. To calculate the susceptibility, we 
imagine evaluating the free energy in a large, finite box 
of linear size L. We differentiate with respect to the 
field, and take i? ^ to get the susceptibility. The limit 
L — > cx) is taken at the end. This gives 



AxiT) 



2nasUt(giiB) 



dk 



^duj 



92 



1 



2tt 1 - e-l^f''^ duj^ cj2 + D2fc4 



^- dk 



^duj 92 



;fc2 



27r duj^ V ef"'''^ - 1 



1 



DU^ 



It is assumed that the fc-integral is cut-off at the lower 
end by the inverse system size L~^ and at the upper end 
by the inverse mean free path. In going to the second 
line, we have performed an integration by parts twice. 
We may now proceed exactly as for the specific heat cor- 
rection above. We first replace the frequency integral by 
one that runs over positive lu alone, and integrate over 
the angular components of the momentum to get 



hut{g^J.BYcrs 



92 



5^2 I ^[jhu _ I 



Dlk^ 



(We denote Jj^ ^ — J^-i dk J^ dw). The fc-integral is 
infra-red convergent, and may be performed first as be- 
fore. The singular contribution may be evaluated exactly 
as for the specific heat above, and yields the result quoted 
in Section ^. 



3. Electron tunneling density of states 

Here we show that the electron tunneling density of 
states is, under certain approximations, essentially the 
same as the d-particle tunneling density of states. For 
concreteness, we consider a lattice electron Hamiltonian. 
The single-particle electron density of sta,tes{Nc{E)) is 
related to the electron Green's function Q by 



1 



Nc{E) = hngu{E + iT]), 



(B2) 



where the overline denotes averaging over the disorder 
and i is a site index on the lattice. The Green's func- 
tion Q may be obtained by analytic continuation from 
imaginary frequencies: 

Guiiuj) = (c,|(w)c,|(a;) + (T^i)). (B3) 

Transfoming to the d-fields, the right-hand side becomes 

{d,^{uj)d,^{uj) - d,i{-uj)d,i{-iuj)) . 

Defining the d-particle Green's function Ga^aaii^) — 
{dia(uj)dia{uj)) , this becomes 



Now spin SU{2) invariance requires Ci|(aj)ci|(a)) = 
Cii{uj)cii{uj) which implies that 

Gjj,tt(icj) = -G^^^ll{-iuJ). (B4) 

We therefore have 

GuiiLu) ^2Gii^ijij{iLu). (B5) 

The d-particlc tunneling density of states is 

Nd{E) = -hm {G^^,^^{E + iri) + Gu.ii{E + ifi)) . 

(B6) 



Using Eqns. B4 and B5, we may rewrite this as 

^Im (gu{E + ii^)-gu{-E-i7j)'^ (B7) 

(B8) 



NdiE) 



-{N,{E)+N,{-~E)). 



If wc now finally assume that asymptotically close to the 
Fermi energy, there is a statistical particle-hole symme- 
try for the electrons, then Nc{E) « Nc{-E) as E ^ 0. 
We then have NdE) -> Nd{E) as E ^ 0. 

The methods of Ref. |0] can be used to show in a 
straightforward way the existence of a VE singularity 
due to interactions in Nd{E) in the superconducting ther- 
mal metal in three dimensions. As we noted in Section 
III , Nd{E) has singularities due to quantum interference 
as well. We therefore have, for the electron density of 
states Nc{E), at low energies 



NciE) - Nc{0) 



E, 



(B9) 



with contributions from both quantum interference and 
interaction effects. 
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